QMUL-PH-07-11 
RH-04-2007 



Type I singularities and the Phantom Menace 



Tapan Naskar a and John Ward b ' c H 



a IUCAA, Post bag 4, Ganeshkhind, Pune 411007, India 

b Center For Research in String theory, Department of Physics, 
Queen Mary University of London, Mile End Road, London, El 4NS, UK. 

c Science Institute, University of Iceland, Taeknigardi, Dunhaga 5 IS- 107 Reykjavik, Iceland. 

Abstract 

We consider the future dynamics of a transient phantom dominated phase of the universe in 
LQC and in the RS braneworld, which both have a non-standard Friedmann equation. We find that 
for a certain class of potentials, the Hubble parameter oscillates with simple harmonic motion in 
the LQC case and therefore avoids any future singularity. For more general potentials we find that 
damping effects eventually lead to the Hubble parameter becoming constant. On the other hand in 
the braneworld case we find that although the type I singularity can be avoided, the scale factor still 
diverges at late times. 



1 Introduction 



Observations suggest that the universe is currently dominated by a dark energy phase, which accounts 
for the apparent acceleration pQ. However this cannot be reconciled within conventional cosmology, and 
therefore provides an important theoretical challenge |2j • One possible explanation for this dark energy 
is that there exists a non-zero cosmological constant, implying that we are currently in a de Sitter phase. 
The origins of such a constant have been the focus of much recent work in flux compactifications of 
superstring and M-theory (see [H] and the references therein.). An alternative proposal is that the energy 
density of the universe is currently dominated by a phantom scalar field, which has a canonical scalar 
Lagrangian albeit with the wrong sign in the kinetic term |10j . This leads to a violation of the weak 
energy condition such that p < p. The physical origins of such a field, however, remain problematic at 
best since the Hamiltonian is unbounded from below which means that a stable vacuum state may not 
exist upon quantization of the theory. In addition is was shown that phantom fields must have minimal, 
or zero, interactions with other matter - and that they themselves must only interact gravitationally [7]. 
However see for a potential resolution to the phantom problem. 



One resolution to this problem is that the phantom fields may only exist temporarily, which is what 
we will assume in this note (although we should point out that these kind of fields can naturally arise 
within certain string theories, most notablytype II* string theory [5]). Let us imagine that the kinetic 
term in the Lagrangian has some moduli space metric Gij (<^>) multiplying the derivative terms such that 
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we have a term Gijd^d^ft . Usually one takes the metric to be positive definite, however let us assume 
that there is some region of moduli space where it has a negative eigenvalue. Thus any fields entering 
this region will have negative kinetic energy, at least for the duration of time that the field find itself in 
this region of moduli space [5] ■ It is this region that we shall consider in this note which therefore implies 
that we expect the current dark energy phase to be a transient phenomenon. See for related issues. 

Of course the scalar field may be remain in this region for a long time, and thefore the accelerating phase 
may appear to be long lived. 

In a typical Friedmann-Robertson- Walker (FRW) model exhibiting a dark energy, or phantom phase, 
such a universe will contain a future singularity. In this type of singularity, referred to as a Type I 
singularity [3], both the energy density p of the dark energy and the Hubble parameter diverge within 
a finite time. Clearly this is undesirable from a physical perspective, and so many authors have tried 
to resolve these future singularities [15] . Although there has been much progress in resolving initial 
singularities in string theory, future singularities arc problematic. This is because the low energy solutions 
of pcrturbative string theory simply match onto the standard FRW cosmologies. Therefore we are 
motivated to move beyond this approach to search for other solutions. Of course we may find additional 
corrections once a fully consistent non-perturbative description has been realized. 

An alternative approach to quantum gravity is via loop quantum gravity. In this theory non- 
perturbative effects lead to p 2 corrections to the standard Friedmann equation and thus allow us the 
possibility of resolving any future singularities [H H31 HH [H] . It is this non-canonical dependence upon 
the energy density that leads to interesting physical behavior, particularly when the energy density is 
near the critical density p c . As such backgrounds of this type deserve further investigation. Recently it 
was observed that these non-perturbative corrections can lead to the future oscillation of H between finite 
values, implying that neither the scale factor nor the energy density diverge. Instead both parameters 
become oscillatory and bounded at finite value |22| . thus avoiding any future singularities. 

Another interesting correction to the Friedman equation can be seen in studies of brancworld cosmol- 
ogy |17| . where there is a p 2 correction which has the opposite sign to the loop background. In transpires 
that there is an interesting and highly non-trivial duality between these solutions |14] suggesting that the 
two models are related to one another. 

In this note we will study the qualitative differences between the loop and brancworld inspired solu- 
tions, and their potential consequences for the evolution of the universe. We will see that in the loop case 
the Hubble parameter behaves as a damped harmonic oscillator with frequency uj, however in the brane 
case the oscillation frequency rapidly becomes complex and the Hubble parameter is frozen at a constant 
value Although this avoids the type I singularity, the scale factor continues to diverge and therefore we 
do not resolve the singularity issue. 

Note: Upon completion of this work we learnt of the paper by [27| which investigates a similar issue. 
The conclusions reached there are similar to those presented here. In both cases we emphasis that the 
Hubble parameter can become oscillatory, although our analysis shows that this will only be true for 
a specific class of phantom potentials. We also highlight the similarity between this model, and that 
coming from a brancworld motivated one. 
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2 Phantom Scalar Fields and Loop Backgrounds. 



The energy density p and pressure p of a phantom field are defined by the following components of the 
energy momentum tensor 

p = -\4> 2 + v 

P = (2-1) 

and the equation of motion 

4> + 3Hj) = V'. (2.2) 

This is enough to specify the dynamics of the minimally coupled phantom field coupled to standard 
Einstein gravity. The standard Friedmann equation, in terms of the reduced Planck mass M 2 = 8irG, is 
simply 

H 2 



however non-perturbative loop quantum gravity effects modify this to rea 

3M 2 V p ( 



tf 2 = ^(l-^) (2-3) 



where V = dV/d(j> and p c = y/3/ (16irj 3 G 2 h) is the critical energy density. Here 7 is the Barbero-Immirizi 
parameter which arises from the Poisson bracket of the conjugate connection c and the triad p, through 
the relation {c,p} = 8ttGj/3. Numerically 7 is typically fixed to be small. We refer the interested reader 
to [12] for a detailed account of the quantization procedure for such a theory. Combining these we can 
write the derivative of H with respect to time 



P + P 



-^777^(l--l (2-4) 



As usual in these problems it is easier to work in terms of dimcnsionlcss variables, so we make the following 
definitions 



/3M 



2 



V 

V=— (2.5) 



3M2 Pc 

and using these dimcnsionlcss quantities we see from (|2.3l I2.4|) that the dimcnsionlcss equations of motion 
become. 

H 2 . + 



2 V dt J I \ 2\dt 

dV_ fl +3 H^ 
dj> dt 2 dt 



3 This equation also arises in braneworld constructions with an additional time-like direction, as shown in Shtanov and 
Sahni El- 
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We want to re-write these expressions in an autonomous form, so it is convenient to define the following 
variables 



dt 



X 



r 

v #' 



y= vv 
vv" 



y/2 



And thus the autonomous equations of motion can be written 

dx 

dt 
dy 

dt 



= -3Hx - Xy 2 
-\\xy 



= a i — r 



dX 
dt 

Additionally we see that the Friedmann equation in dimensionless variables becomes simply 



(2.7) 



(2.8) 



H 2 = \-\*? + y 2 )[l + \x 2 -y 2 



(2.9) 



The non-negativity of this last expression implies that (x 2 /2 — y 2 ) G [—1,0] or in terms of the variable 
x 2 we see that this must lie in the closed set x 2 G [2(y 2 — 1), 2y 2 ]. The variables are constrained, and not 
independent which is why the Hubble parameter never becomes imaginary even when x = 0. This can 
be demonstrated graphically using a Monte Carlo simulation and shown in the figure below 
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Figure 1: Phase space scatter plot of the {x,y} plane generated using a Monte Carlo algorithm. The 
boundary (blue lines) defines the region of allowed points. 
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As we move away from the origin we see that x becomes proportional to y as expected from our 
analytic argument. Since y 2 is the dimcnsionlcss potential for the phantom field this indicates that we 
should find scaling solutions when y 2 >> 1, so that we can neglect all 1/y 2 terms as being negligibly 
small. 

Importantly we can write the acceleration of the Hubble parameter in the following form using (|2.6[) 

*» - -4 +2 .-w)*-3v(* + .-w) 

= -u?H-^ (2.10) 

Where we have defined 

lu 2 = 9x 2 ^1 + 2x 2 - 2y 2 ^j (2.11) 
v 2 

7 = 2A— (2.12) 

x 

We see that H is a damped, harmonic oscillatory function with frequency uj and damping factor 7. This 
is an important result, since this implies that H can oscillate and therefore avoid future singularities. The 
damping factor scale is set by A, and decreases as y increases, therefore one would expect the damping to 
die away as a function of time. The frequency is an ever increasing function of y provided that y 2 > 3/2, 
which is a reality condition. One can see that the frequency will never become complex, because the 
value of y is constrained by the value of x. Again we can confirm the validity of said assumption using 
Monte Carlo techniques 



Therefore it would appear that at late times H oscillates more rapidly, becoming a simple harmonic 
oscillator. Schematically we see that the Friedmann equation looks like H 2 ~ p(l — p), and so we deduce 
that p must also be an oscillatory function. 

We can define an effective equation of state for the phantom field in the usual manner 




(2.13) 



which essentially tells us that phantom dominated solutions require x 2 > since all dependence on the 
potential vanishes when we look at violations of the weak energy principle. Despite this, we know that 
the kinematics of the phantom field are determined by the form of the scalar potential. Therefore we 
should try to analyze the equations of motion with this in mind. There are essentially three cases of 
interest which we label as i) A = const, n) A — ► asymptotically and Hi) A — > 00 asymptotically which 
we will discuss in turn. 



2.1 A = const. 



This case implies three the existence of three different possibilities. Firstly we could take A 2 = 0, which 
essentially implies that V = or that the potential is a constant. The second possibility is that x = 0, but 
this means that the phantom field is constant and doesn't permit dynamical solutions. The third and final 
possibility is that T = 1 which uniquely fixes the form of the potential to be of the form V((f>) = A cxp(a<fr) , 
where A, a are positive or negative constants of the appropriate dimension. A nice illustration of the 
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Figure 2: Plot of the square of the frequency as a function of y, using the causality constraint. The 
scattered points represent the allowed values of uj 2 . It is noticeable that there is a band of physical 
solution which increases in width as y also increases. 

expected oscillatory behavior is shown in figure 3, where we use the potential V ~ e - ^ to illustrate the 
general behaviour. Note that the Hubble parameter initially undergoes a bounce, before settling into an 
oscillatory phase 0. The amplitude of the oscillations varies rapidly between — 1/2 < H < 1/2. Similarly 
the energy density settles into an oscillatory phase, after initially increasing towards its maximal value. 
It is important to note that the Hubble parameter is approximately constant or slowly decreasing for 
much of the initial evolution, therefore there will be virtually no corresponding change in the amount of 
gravitational wave production. Eventually the Hubble parameter reaches zero, at which point is moves 
into its oscillatory period (with a fixed amplitude), but oscillating about the zero point. This will therefore 
lead to an oscillatory signal for the gravitational waves produced, with an ever decreasing period. This 
oscillation is expected to be eternal [57]. Physically the bounce is due to the way the energy density is 
distributed into kinetic and potential energies. The potential energy of the field is always positive, and 
initially dominates the negative kinetic energy of the rolling field, which is why the energy density is an 
initially increasing function. However eventually the field starts to roll rapidly and therefore the kinetic 
energy contribution starts to become dominant. Oscillation in the velocity thus makes the total energy 
density oscillatory, which in turn forces the Hubble parameter to oscillate about the bounce point. 

Let us initially assume that x c = 0. In fact the only fixed point solution can be see to be x c = 0, y c = 
- which is essentially the most trivial case because the energy density is zero. We can actually solve the 
equations of motion in this instance. We see that y is constant, and that the boundary term for the 
velocity must satisfy the equation xq = Xy 2 {t — to). However the equation of state at this point must be 
evaluated as a limit Lim Xiy= of20 in order to avoid singularities. In fact this particular limit can be shown 
to be non-existent when one applies the tests for multivariable functions, therefore the equation of state 

4 This was also discussed in |25| . 
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Figure 3: Oscillation of the Hubble parameter and energy density when we take A = const, using a 
potential given by V — e~^. 



is undefined in this instance. 



Let us now consider solutions with y c = 0, which implies that the fixed point is again at (0, 0) if we 
assume that the phantom field doesn't have complex velocity. The more interesting case is to set A = 0, 
which implies that the following critical points exist (0, y), (±\/2y, y) and (±-»/2(y 2 — 1), y), with y taken 
to be arbitrary. The point (0,y) admits a constant equation of state - namely f2^ = —1 consistent with 
a dark energy dominated phase rather than a phantom phase. The point (±^/2y, y) leads to a divergent 
equation of state with Q$ — > — oo for all non-zero positive y. The final fixed point actually leads to the 
solution fl,p = 1 — 2y 2 and therefore to a whole range of phantom like behavior provided that y 2 > 1. 
Solutions to the equation of motion can be obtained numerically, however they show that x — > as a 
function of time, and that — > — 1 from below. 



Let us now examine the stability of the solution by considering perturbations around the critical 
points of the dynamical equations, where we assume the perturbation is small 

x = x c + 6x, y = y c + Sy. (2-14) 

We wish to determine the perturbed equations of motion to leading order. The only non-trivial expression 
we need to determine is the perturbed Friedmann equation. If we denote H 2 = H 2 (x c ,y c ), i.e the 
Fricdmann equation evaluated at the fixed points - then we can write the perturbed solution as follows 

H~±H e (l + ^ fa - 2 ^ffi 2 - 1 -^ + ...) (2.15) 

We now insert this into the original equations of motion and expand to linear order in the perturbations. 
We can write the resultant differential equations for the perturbed solutions as a matrix equation where 
we define a perturbation vector 5Z = (5x, dy), and a perturbation matrix M satisfying 

|(<5Z)=M(«5Z), (2.16) 
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whose eigenvalues determine the stability of the solution in question. A brief calculation yields the 
following solution for the eigenvalues. Let us first define the function F(x, y) through the relation 

F{x c ,y c ) = 9Ht-12H 3 c \y c ~12H c \x c y c {x c + 2y c )[2y 2 c ~l~xl}+ (2.17) 
+ 9[2y 2 - 1 - x\\ 2 x\(x c + 2y 2 ) 2 + H 2 {l8x c (x c + 2y c [2y 2 c - 1 - x 2 ] + X 2 y c (x c + Ay c ))} 

in terms of which we can write the eigenvalues as follows 

/?± = -3H 2 - 3x c (x c - 2y c )[2y 2 - 1 - x 2 c ] - 2H c Xy c ± F x l 2 {x c , y c ) (2.18) 

The system can be regarded as being stable (at least to linear order in the perturbations) when both 
eigenvalues are negative definite In the first instance where we have (x c ,y c ) = (0,0) the eigenvalues are 
in fact undefined, since there is no possible path through phase space that leads to a convergent solution. 
This should have been anticipated from the form of £1$, which cannot be defined in this limit. 

The non-trivial case where we have A = leads to different solutions. For the point (0, y c ) we 
find that the eigenvalues are = 0, —6H 2 implying that this point is a stable node provided we take 
the positive sign for the Hubble parameter. The remaining fixed points correspond to (±y/2y c , y c ) and 
(±y^2(y| — 1), y c ) and yield zero eigenvalues. This means that they are both saddle points in phase 
space. We summarise the results in this section in the following table: 



Validity of A 


Fixed points 


Eigenvalues 




A = arbitrary 
A = 
A = 
A = 


(0,0) 
(0,2/) 
(±V2y,y) 
(±V2(y 2 -l),y) 


undefined 
0± = 0, -6H 2 

P± = o,o 

= 0,0 


undefined 

^ = -1 
fltf, —> — oo 
% = 1 - 2y 2 



Table 1. Summary of the fixed points and equation of state for the case of constant A. 

We also plot the phase space trajectories (Figure 4) for x and y in the specific case where we select 
A = 0.5, where we have been careful to ensure that the trajectories always remain within their regime of 
validity. This is representative of almost all the small, but constant A trajectories. 

What is clearly visible from this plot is the trajectories appear to converge at the point y ~ 1.05. 
Note that the top line in the plot corresponds to the curve y = yl + x 2 /2. The other curves tend to 
approach this one at late times, despite having vastly different initial conditions. Of course this will only 
be an approximate relation since A is small, but non-zero. If we insist on setting this to zero, then we see 
that the trajectory is simply a straight line which wc have opted not to plot. A similar thing happens 
when we take the converse limit and assume A >> 1. The resulting trajectories simply collapse onto a 
straight line as you increase A to larger values. 



2.2 A^0. 

There are a whole class of potentials that satisfy this limit. The simplest being an inverse potential 
of the form V((j)) = Vo</>~", where Vq is an appropriately normalized mass parameter and the phantom 
is an increasing function of time. Other examples along these lines correspond to V ~ Vbe 1 /^ and 
V ~ Voe M $ I 2 . In fact the A = constant solutions discussed in the previous section can be thought of 
as being instantaneous snapshots of these solutions. Physically we expect the phase space trajectories, 
such as those displayed in figure 4, will tend to collapse onto a single curve as this limit is approached. 
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The equations of motion in this approximation, provided we ensure that T ^ 1, become 

dx ~ dy dX , 

Clearly from the first of these equations we see that x < when x > 0, and also that x > when 
x < which suggests that the attractor solution is x — » asymptotically. The other expressions again 
imply that y and A both tend to constants at late times. How does this affect the expected oscillatory 
solutions? We see that the frequency will tend to zero at late times regardless of the specific value of 
A and therefore there will be no oscillation of the Hubble parameter. One can see this intuitively since 
the damping factor diverges in this limit (because x eventually settles at the point x = 0), provided that 
y =/= 0, implying that the solution is overdamped. Note that A is tending to zero, but may remain small 
but non-zero since it must tend to a constant at late times. 



2.3 A — > oo 



This condition is much harder to investigate because of the equation of motion for A. In order for this 
to be an increasing function of time we must ensure that either T < 1 with x > 0, or that T > 1 with 
x < 0. The first of these cases is not consistent since one can sec that this condition implies that A < — 1, 
therefore only the second solution is physical. In fact this implies that the potential must follow a power 
law solution. From the equations of motion we see that the damping factor increases with time, eventually 
killing the oscillatory behavior of the solution. 
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A nice example of potentials satisfying this condition are power law solutions where V = M^~ n <\> n . 
This leads to increasing T when n > 1. Figure 5 illustrates the oscillation of the Hubble parameter 
for the </> 4 case, where we work explicitly in Planck units. Again we see the bounce solution for the 
Hubble parameter - and also the initially increasing energy density of the phantom field, before both 
become oscillatory. The initial increase of both parameters is relatively quick when one compares this 
to the solution plotted in figure 3, and there are significantly more oscillations, each with a reduced 
period. This is therefore a falsifiable signal of the A — > oo behaviour. Interestingly the Hubble parameter 
decreases fairly steadily at early times indicating that the universe continues to expand, but with a 
decreasing velocity. This should be contrasted with the plateau phase of the A ^constant model. 




-0.6 I 1 1 1 1 1 1 1 1 
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Figure 5: Oscillation of the Hubble parameter and energy density assuming a simple <fi 4 potential. Note 
that the initial increasing phase is relatively short lived, giving way to many oscillations with an ever 
decreasing period of oscillation. This suggests significantly different behaviour than the A = const, 
solution illustrated earlier. 



3 Phantom Braneworlds. 



Recently |14j it has emerged that there is a relationship between loop quantum gravity effects and the 
Randall-Sundrum (RS) (16] brancworld. Recall that in the braneworld scenario there is a branc localized 
in a 5d AdS bulk spacetime, where the extra dimension is taken to be large. This was shown to naturally 
confine Einstein gravity to the world- volume of the brane [17] . Although motivated by compactifications 
of the strongly coupled heterotic string |20j , the braneworld scenario is not derivable from a fundamental 
theory. However it does provide us with an important testing ground for model building, and in the light 
of the duality with LQC may offer new insight into the problems of loop quantum gravity. It is not the 
purpose of this note to review the braneworld scenario, we refer the interested reader to [18] . 

In fact - the RS Braneworld has the same autonomous equations of motion as in the loop case |17j , 
except for the fact that p c is given by a slightly different expression involving the tension of the brane. 
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Our analysis therefore requires us to consider the branc tension as contributing to the critical density. 
There is also a change of sign in the second term in the Fricdmann equation so that schematically we 
have H 2 ~ p(l + p), and thus the reality bound on x becomes x 6 ±^/2[y. \J\ + y 2 ]. 




The expression for the acceleration of the Hubble parameter remains the same except that we must now 
define the frequency as follows 

uj 2 -> 9x 2 (1 + 2y 2 - 2x 2 ) (3.2) 

This has important consequences for the behavior of the solution, as we can see that the frequency is 
only R in a finite domain, and that it becomes complex at late times. This means that H will initially 
try to oscillate, but as time evolves it reaches the critical point and freezes out at some constant value 
0. This is completely different to the loop behavior and offers us some falsifiable prediction. One can 
see this since the condition for uj to be complex is when 2x 2 > 1 + 2y 2 , which can lie within the bound 
coming from the Fricdmann equation, where we must ensure x 2 < 2(1 + y 2 ). Therefore there is a region 
of phase space where H 2 is real, but the frequency of the oscillation becomes complex. 



We can analyze the attractor solutions in the same way as we did for the loop case (see also |19j). 
Since the formalism was already discussed earlier we will be brief in our description. Let us first consider 
the case of constant A as before. The equations of motion admit four critical solutions once again, which 
are given by (0,0) for arbitrary A and (0,y c ), (±V2y c ,y c ), (±^2(1 + y 2 ),y c ) when we take A = 0. Even 
without analyzing the stability of such solutions we can use the effective equation of state to understand 
the physics. The first solution again leads to an undefined equation, for the same reason as before. The 
second solution, with y corresponding to a flat direction, yields f2^, = — 1 which is pure dark energy, and 
doesn't represent a phantom solution. The final solution yields an equation of state given by fl^ = l + 2y 2 
which is always greater than zero. In fact it is only the third of the above solutions leads to phantom 
type behavior. 



To analyze the stability of the solutions we must again calculate the eigenvalues of the perturbation 
matrix. Once again we define the function F(x c ,y c ), which this time becomes 

F{x c ,y c ) = 36H*-48H?\y c + 24[l + 2y 2 c -x 2 c ]H c \y c {x c + y c ) + (3.3) 
+ 9[1 + 2y 2 - x 2 c ] 2 (x c + y c f - 4H 2 {9[1 + 2y 2 c - x 2 c ](x c + y c ) - \ 2 y c (x c + 4y c )} 

and therefore we can write the corresponding eigenvalues as 

a 3(z c - y c )[2y 2 - x 2 } + H C (-6H C - AXy c ± F 1 / 2 {x c , y c )/H c ) 

,3 ± = . (3.4) 

The point (0, 0) is the same as in the loop case, therefore we will not comment on it further. The second 
fixed point is again at (0,j/ c ), however this time the analysis is complicated by the form of the Hubble 
parameter. We see that the eigenvalues become 

/?=-!- (-2H 2 C -y c - 2yl ± \2H 2 -y B - 2y 3 c \) . (3.5) 

Now if we take H c to be positive definite then we find that the two eigenvalues are both negative functions 
of y c , and therefore correspond to stable points in the phase space. Conversely if we chose the minus sign 
for the Hubble parameter, then we see that both eigenvalues are positive definite, and this corresponds 
to a repulsive fixed point. The remaining two fixed points correspond to the zeros of the Friedmann 
equation, and much like the loop case lead to eigenvalues that are identically zero. This indicates that 
they are saddle points in phase space. We again summarise these results in the following table 



3 The oscillation continues but in the complex plane 
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Validity of A 


Fixed points 


Eigenvalues 




A = arbitrary 
A = 
A = 
A = 


(0,0) 

(O.tf) 
(±V2y,y) 

(±v/2(y2_ 


undefined 

0± = -2H 2 - 2/(1 + 2y 2 ) ± \2H 2 - y(l + 2y 2 )\ 

P± = o,o 

= 0,0 


undefined 

« = -1 
— ► — oo 
ft = 1 + 2y 2 



Table 2. Summary of the fixed points and equation of state for the case of constant A for the 

braneworld case. 



4 Conclusion 



In this note we studied the dynamics of a phantom driven universe both in a loop and braneworld 
background. We found the Hubble parameter H and the dark energy density p remain finite in future 
and hence avoid a type I singularity for the loop case, however in the braneworld case we see that the 
scale factor can still diverge due to the fact that the frequency of the oscillation becomes complex. We 
analyzed the fixed points of the equations of motion in both cases and found that there were two non-trivial 
solutions, and two trivial (in the sense that they correspond to H = 0) solutions. The first non-trivial 
solution is undefined in both cases, whilst the second is a stable fixed point with the appropriate choice 
of sign. The trivial solutions both correspond to flat directions, even though that they are located at 
different points in the (x, y) phase plane. 

Our study of the loop background was for a general phantom potential, however we saw that the 
equations of motion impose restrictions on the allowed form of such potentials. The most interesting 
solution corresponds to the A = const, of which the A = solution is a subset. Potentials satisfying this 
condition lead to simple harmonic motion for H in the far future, and are therefore capable of resolving 
all future singularities. On the other hand, the asymptotic cases where A — > 0, oo should appear to be 
oscillatory only for a finite time. Eventually the equations of motion force the damping factor to dominate 
at late times, and the system will be fixed at some value. Of course this bounce appears to pass through 
the point H = which is where we would anticipate quantum gravity effects to become important. 

In the braneworld case we can again use A to parameterise the solutions, however the Hubble parameter 
quickly freezes out to become constant because of the reality constraint on the oscillation frequency - 
even when A = 0. This means that this the braneworld case is similar to the late time behavior of 
the asymptotic loop case. In the asymptotic cases A — ► 0, oo we see that again the damping factor will 
dominate at late times, and therefore kill off all oscillatory solutions, exactly as in the loop case. Therefore 
we conclude that only when we have constant A can we be sure of avoiding all future singularities. In the 
other cases it is likely that the type II singularity (a — > a s , t — > t s , p — > p s ) could still occur [3J, although 
the sudden singularity [3] is regarded as being a far weaker singularity. 

The rapid oscillation of H and p in the loop case leads to some interesting cosmological phenomenon. 
This may cause late time particle production and therefore provides a useful falsifiable prediction com- 
pared to the braneworld scenario. There will certainly be an observable gravitational wave signature due 
to the oscillatory nature of H, however the precise physical predictions are dependent upon the choice 
of the potential, and therefore depend heavily upon the physical origin of such a phantom field. In both 
cases we have assumed that the evolution of the universe is governed by a transient phantom phase, and 
that eventually this phase will come to an end. However this may be problematic, since it was argued in 
[26] that one cannot find trajectories where the eigenvalues of the metric tensor take the opposite sign H. 

6 Many thanks to A. Vikman for pointing this out. 
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One way to avoid this is to consider a large number of phantom fields, in which case it may be possible 
for the moduli space metric to acquire an opposite sign. A detailed analysis of this conjecture would 
certainly be most welcome. However it is important to understand the differences between the brane and 
loop solutions during this temporary phase of phantom domination. 

Assuming that the phantom phase is transient allows us to consider more general corrections to 
the Friedmann equation, such as those arising in Cardassian models (53] or the DGP model [21] ■ It is 
important to understand the distinguish between the different physical predictions of these theories in 
order to understand whether all future singularities can be resolved. 
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